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ABSTRACT 

We consider dark masses measured from kinematic tracers at discrete radii in galaxies for which baryonic 
contributions to overall potentials are either subtracted or negligible. Recent work indicates that rotation curves 
due to dark matter (DM) halos at intermediate radii in spiral galaxies are remarkably similar, with a mean ro- 
tation curve given by log 10 [V c .DM/(kms" 1 )] = 1.47^[g + 0.51og 10 [>/kpc]. Independent studies show that while 
estimates of the dark mass of a given dwarf spheroidal (dSph) galaxy are robust only near the half-light radius, 
data from the Milky Way's (MW's) dSph satellites are consistent with a narrow range of mass profiles. Here 
we combine published constraints on the dark halo masses of spirals and dSphs and include available measure- 
ments of low surface brightness galaxies for additional comparison. We find that most measured MW dSphs lie 
on the extrapolation of the mean rotation curve due to DM in spirals. The union of MW-dSph and spiral data 
appears to follow a mass-radius relation of the form Muu{r)/ Mq = 200+ 2 2o(r/pc) 2 , or equivalently a constant 
acceleration #dm = 3 +f x 10~ 9 cms~ 2 , spanning 0.02 < r < 75 kpc. Evaluation at specific radii immediately 
generates two results from the recent literature: a common mass for MW dSphs at fixed radius and a con- 
stant DM central surface density for galaxies ranging from MW dSphs to spirals. However, recent kinematic 
measurements indicate that M31's dSph satellites are systematically less massive than MW dSphs of similar 
size. Such deviations from what is otherwise a surprisingly uniform halo relation presumably hold clues to 
individual formation and evolutionary histories. 

Subject headings: galaxies: kinematics and dynamics — galaxies: spiral — galaxies: dwarf — galaxies: fun- 
damental parameters — (cosmology:) dark matter 



1. INTRODUCTION 

Considered within the framework of Newtonian mechanics, 
the observed motions of stars and gas imply that the lumi- 
nous components of galaxies are embedded within extended 
halos of unseen material. Absent the ability to measure in- 
visible masses directly, for a given galaxy one measures dy- 
namical and luminous masses independently and then iden- 
tifies dark matter simply as the difference between the two: 
Mum = Md yn -M\ um . Here we shall investigate the properties 
of dark matter halos only insofar as they are accessible via 
the direct application of this equation, without adopting any 
particular halo model. 

This task is made difficult by the fact that there are no galax- 
ies for which both Md yn and M\ um are easily measured. While 
the circular motions of stars and gas in spiral galaxies provide 
a relatively clean measure of the dynamical mass profile, sub- 
traction of the luminous contribution suffers from systematic 
uncertainties, primarily those related to stellar mass-to-light 
ratios. One sidesteps this problem, or at least exchanges it 
for others, by considering dwarf spheroidal (dSph) galaxies, 
for which baryonic components contribute negligibly to inter- 
nal gravitational potentials. However, because dSphs are sup- 
ported by stellar velocity dispersions instead of rotation, esti- 
mation of their individual mass profiles requires large kine- 
matic data sets and typically em ploys strong modeling as- 
sumptions (see discussion by, e.g.. lPryor & Kormendyl[l990l : 
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IWilkinson et alj|2002t iGilmore et al"1l2007l) . Larger elliptical 
galaxies (see iNapolitano et al.ll2010l ) combine the most chal- 
lenging aspects of both regimes, as their significant stellar 
masses are supported in large part by random motions. 

Nevertheless, recent work indicates that dark masses are 
constrained reasonabl y well at intermediate r adii of both spi- 
ral and dSph galaxies. iMcGaugh et all J20071 "M07") use ro- 
tation curve data for a sample of 60 spirals and discard data 
points interior to r < 1 kpc, where dynamical masses can be 
affected by noncircular motions. After subtracting baryonic 
masses, M07 find that the rotation curves due exclusively to 
dark matter halos lie nearly on top of each other, with mean 
rotation curve 

logiotVWCkms- 1 )] = 1.47tli + 0.51og 10 [r/kpc]. (1) 

Working at smaller scales. IWalker et al.1 (120091 "W09") use 
kinematic data for eight Milky Way (MW) dSphs to show 
that for a wide range of plausible halo shapes and velocity 
anisotropies, the estimated mass within the projected half- 
light radius is approximately (subject to validity of the as- 
sumptions of spherical symmetry, dynamical equilibrium, flat 
velocity dispersion profiles (v 2 ) 1//2 (/?) = er, and negligible con- 
tributions to measured velocity dispersions from unresolved 
binary motions) 
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Subsequently lWolf et alJ (120091) have provided an analytic ar- 
gument for why such an estimate is insensitive to anisotropy. 
Since dSph kinema tics tend to b e dominated by dark matter 
at all ra dii (Section | 2.3l see also Mateo 1998 and references 
therein: IWalker et alj|2007h . Equation [2] provides a crude es- 
timate of the dark mass within r^if for any dSph for which 
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measurements of and a are available. W09 apply Equa- 
tion|2]to data for 28 classical dSphs and ultrafaint satellites in 
the Local Group and find that, allowing for scatter by a factor 
of less than two in normalization, the ensemble of dSph data 
is consistent with a mass profile of the form M oc r lA±0A . 

Thus M07 and W09 independently note the apparent self- 
similarity of spiral and dSph dark matter halo profiles, respec- 
tively. Here we combine data from those studies in order to 
examine whether the mean rotation curve attributable to dark 
matter in spiral galaxies extends to the small radii character- 
istic of dSphs. For fu rther comparison we a lso in clude re- 
cent measurements by iKuzio de Naray et alJ (120081) of rota- 
tion curves in low surface brightness (LSB) galaxies. The 
combined data set samples dark matter halos over the range 
0.2 < r < 75 kpc. 

2. DATA 
2.1. Spirals 

T he M07 data consist of HI rotation curves compiled 
by ISanders & M cGaugh (120021 observational references 
therein ), trimmed for qu ality to a formal accuracy of 5% or 
better ( McGaughl 120051) . and restricted to radii larger than 
r > 1 kpc because systematic uncertainties (e.g., noncircular 
motions, beam smearing) dominate the deconvolution of ro- 
tation curves at smaller radii. The final sample contains 686 
independent, resolved rotation velocity measurements for 60 
galaxies, spanning radii 1 < r < 75 kpc. This sample covers 
virtually the entire range of spiral properties, including circu- 
lar velocities 50 < Vf < 300 km s , scale radii 0.5 < rj < 13 
kpc, baryonic masses 3 x 10 s < Mj < 4 x 10 n M Q , central 
surface brightnesses 19.6 < /io b < 24.2 mag arcsec 2 , and gas 
fractions 0.07 <f g < 0.95. 

In order to isolate the contribution of dark matter to an ob- 
served rotation curve, M07 subtract contributions from stel- 
lar (bulge, disk) an d HI components. For this task M07 
use GIPSY software (Ivan der Hulst et aO 19921) to convert nu- 
merically the observed luminous distributions into rotation 
curves. The most critical step in this procedure is the esti- 
mation of the stellar mass-to-light ratio, T* = M*/L*. For 
each spiral, M07 adopt the value of T* obtained by fitting the 
rotation curve with modified Newtonian dy namics (MoND; 
lMilgroml[T98H ISanders & McGaughl 120021) . MoND gener- 
ally provides good fits to the rotation curves of spirals with T* 
as the only free parameter, and the resulting estimates of T* 
sta nd in excellent agreement with popula tion synthesis mod- 
els dBell & de Jongll200UlMcGaughll2004l) . Note that despite 
this use of MoND as an empirically justified fitting formula 
with which to estimate T*, the ensuing dynamical analysis of 
M07 is purely Newtonian. Alternative prescriptions for esti- 
mating T* (e.g., maximum -disk or the population-synthesis 
models of Bell et alJ 120031) give similar results for the dark 
matter rotation curves, albeit with more scatter (M07). 

2.2. Low Surface Brightness Galaxies 

IKuzio de Narav et alJ d2006l [20081 l2QQ9h report high reso- 
lution observations of LSB disk galaxies. These data are two- 
dimensional velocity fields of the inner regions of LSBs, aug- 
mented by long-slit and, at larger radii, HI data. Here we in- 
clude data for the nine LSBs that have m ass models presented 
in the second part of Table 3 from IKuzio de Narav et alJ 
(2008). LSBs are dark-matter dominated down to small radii, 
but their stellar components are not completely negligible 
in the inner kpc. For consistency with the spiral data, we 



adopt the MoND estimates of T* from lde Blok & McG augh 
(1998) where available. These are not available in the cases 
of NGC 4395, UGC 4325, and DDO 6 4, for which we adopt 
popula tion synthesis estimates given by lKuzio de Naray et al. 
(2008). The MoND estimates of are consisten t with those 
of population synthesis models dMcGau gh 2004) with most 
of the velocity attributed to dark matter in either case, so it 
makes little difference which is adopted. The observed veloc- 
ities have also been corrected for pressure support (typically 
~ 8 kms" 1 ); this correction is barely perceptible in most cases. 
We relax the accuracy requirement to 30% in velocity for the 
LSBs in order to retain a reasonable amount of data. 

2.3. Dwarf Spheroidals 

We adopt the dSph data listed in Table 1 of W09 (refer- 
ences to original studies therein; see erratum for updated val- 
ues of rhaif)- These data include projected half-light radii, 
global velocity dispersions and luminosities measured for 28 
objects in the Local Group, including satellites of the MW 
and M31, and two Local Group dSphs — Cetus and Tucana — 
that are unassociated with any obvious parent system. To 
these data we add measurements for eight (six previously un- 
measure d) M3 1 satellites, ad opting velocity dispersions mea- 
sured bv lKaUraietai] d2009h for Andl, Andlfl Andlll, And- 
VI I, AndX and AndX IV, and velocity dispersions measured 
by lCollins et al] J2009I) for And 1X0 and AndXII. 

For each of the 34 objects in the final sample, we apply 
Equation|2]to estimate the dynamical mass, M(rh a if), enclosed 
within the sphere having radius equal to the projected half- 
light radius. We adopt for all dSphs a crude st ellar mass- 
to-lig ht ratio of T* = 1.5 ± l.0M & /L V)Q (e.g., iMartin et all 
120081) . The dark mass is theifl 

M DM (r ha i f ) = M(r half ) - T*^. (3) 

For comparison with the dark masses in spirals we explic- 
itly consider for dSphs the values of MoM('"haif) obtained from 
Equation [3] but because most dSphs are dominated by dark 
matter even at their centers, these values generally are indis- 
tinguishable from the dynamical masses M(rh a it) estimated by 
W09. 

We note that while the velocity dispersions of most 
dSphs are measured from luminous red giants, the published 
dispersions of five ultrafaint MW satellites — Bool, Coma, 
Segue- 1, UMall and Willman-1 — are measured primarily 
from stars near the main-sequence turno f fs of these object s 
(ISimon & Gehal 120071: IMartin et al.ll2007l IGeha etal.ll2009h . 
Since these stars have not yet expanded to become red giants, 
they can still exist as members of unresolved binary systems 
with separations and speeds that are smaller and faster, re- 
spectively, than those available to systems involving red gi- 
ants. Because errorbars on the velocity dispersions of these 
objects do not include systematic errors due to the unknown 
level of contamination by binary orbital motions, masses cur- 
rently estimated for the faintest satellites should be considered 

7 F or Andll, the new measurement of a = 7.3 ± 0.8kms - 1 IKalirai et all 
15003) replaces the value cr = 9.3 ± 2.7kms-' I Cote et al. 1999) listed in Table 
lofW09. 

8 For AndIX, the new measurement of cr = 3.0~^'q km s" 1 'Collins et alJ 
120091) replaces the value cr = 6.8 ± 2.5kms~' IChapman et al.f 2005) listed in 
Table 1 ofW09. 

9 For the Plummer surface brightness profiles adopted by W09, division of 
the total luminosity by 2 3 / 2 gives the luminosity inside the sphere with radius 
equal to the projected half-light radius r^. 



Dark Matter Halos of Spirals, LSBs and dSphs 



3 



as upp er limits. See Olsze wski et all (ll996l) : lHargreaves et alj 
( 1996); Mi nor et all (I20T0F for relevant discussions of the ef- 
fects of binaries on estimates of dSph velocity dispersions. 

3. COMPARING THE DARK MATTER HALOS OF SPIRALS. LSBS 
AND DSPHS 



4 spirals 
. MW dSphs 
□ M31 dSphs 

, Cetus, Tucana 
. LSBs 




FIG. 1. — Top: Circular velocities due to dark matter halos at intermediate radii o f 
spirals (data from Sanders & McGauah 2002: McGauah 2005: McGauah et al. 2007), 
LSBs (data from Kuzio de Naray et al. 2006, 2008, 2009) and within the projected half- 
light ra dii of dwarf sphero i dal galaxies (data f rom Walker et al. 2009 and references 
therein; Kalirai et al. 2009; Collins et al. 2009). Downward-pointing arrows indicate 
five ultrafaint dSphs (Bool, Coma, Segl, UMall and Willi) with velocity dispersions 
measured primarily from faint turnoff stars, for which unresolved binary orbital veloc- 
ities may be significant. Middle and bottom panels give the same information in terms 
of enclosed masses and accelerations, where #dm = G Mdm/V 2 = V^dmA- In all pan- 



els, the solid line indicates the fit by McGaugh et al. (2007) to the spiral data alone: 
log ]0 [V' t . : D]vl/(kms _1 )] = 1.47^5^9+0.5 log 10 [r/kpc], corresponding to a mass-radius re- 
lation of the form M DM (r)/MQ = 200! 2 ™(r/pc) 2 , or g D M = 0.9^ 5 km 2 s~ 2 pc~' (dotted 
lines indicate the formal scatter). Open circles identify the two most extreme outliers — 
Hercules (rhaif ~ 330 pc) and Sagittarius (rhair ~ 1550 pc) — among MW dSph satel- 
lites. Black aiTows indicate the magnitude and direction that a generic dSph would be 
displaced due to the tidal stripping of 99% of the original stellar mass, from N-body 



simulations by Penarrubia et 



"lpping of 
m 120081) 



Figure Q] displays the union of spiral, LSB and dSph data. 
The top panel in Figure Q] plots the rotation velocity due to 
the dark matter halo against the radius where this quantity is 



evaluated. We re-emphasize that these estimates are derived 
independently of any particular (i.e., cored or cusped) halo 
model. Data points for spirals are the same as those in Figure 
3 of M07. For the dSphs we translate estimates of MoMf/haif) 
into circular velocities using V 2 DM = GMum/^- The middle 
panel contains the same information cast in terms of enclosed 
dark mass. Here we plot dSph masses directly and convert 
circular velocities for spirals and LSBs into enclosed masses. 
In both panels, a solid line indicates the fit provided by M07 to 
the mean rotation curve of the spiral data alone, which is given 
by EquationQ]or equivalently by the mass-radius relation 



M DM (r) 



= 200: 



+200 
120 



(4) 



The bottom panel of Figure Q] again conveys the same infor- 
mation, now cast in terms of the acceleration due to the dark 
matter halo, #dm = ^dmA- Since the M07 relation has the 
form M oc r 2 , the corresponding acceleration profile is con- 
stant: 



gDM 



: 0.9!^kmV 2 pc" 



: 3 + l x 10" 9 cms- 2 . 



(5) 



In order to evaluate whether the dSph and LSB data are 
consistent with the extrapolation of the M07 relation at small 
radii, we calculate 



,2 _ 



E 



(log 1 o[VDM,,/(kms- 1 )]-0.51og 10 [r,/kpc]-1.47) 2 

2 2 
CT log 10 [V DM .,/(kms-i)] + CT M07 



(6) 

where residuals are scaled by measurement errors convolved 
with the intrinsic scatter in the M07 relation. The value of 
o"M07 ~ 0.15 reported by M07 reflects the non-Gaussian dis- 
tribution of spiral velocities at a given radius. Here we adopt a 
smaller value of <7mo7 = 0.08 so that the central 68% of spiral 
data points lie within ±omo7 °f tne M07 relation. 

For the 23 MW dSphs we obtain y^/N = 2.2. At present 
the most extreme outlier is Hercules, which has a small ve- 
locity dispersion — and hence a small inferred mass — r elative 
to other MW dSphs of similar size dAden et alj 120091) . At 
least some of the scatter within the dSph satellite popula- 
tions likely results from tidal interactions with the MW and 
M31. Black arrows in Figure Q] ind i cate, according to N- 
body simulations by Penarrubia et al. (2008), the magnitude 
and direction of the track followed by a generic satellite as it 
loses 99% of its original stellar mass to tidal disruption. Tidal 
evolution tends to carry individual objects toward smaller in- 
ferred masses as velocity dispersion decreases more quickly 
than half-light radius. Thus tides may explain the relatively 
small mass estimated for Hercules, whose elongated morphol- 
ogy (axis ratio ~ 3 : 1) may indicate ongoing t i dal disrup- 
tion JBelokurov et alj|2007t [Coleman et al.ll2Q07t ISand et alj 
2009). Lending support to this explanation is the fact that 
the next most severe outlier among MW satellites is the main 
body of Sagittarius, f rom which tidal arm s extend nearly twice 
around the sky (e.g., Maiewski et al. 2003J). If we discard Her- 
cules and Sagittarius as systems that are strongly affected by 
tides, the 21 remaining MW dSphs are reasonably consistent 
with the M07 extrapolation, having y^/N = 1.2. If we fur- 
ther exclude Bool, Com, Segl, UMall and Will, we obtain 
X 2 /N=1.5. 

The 182 points contributed by the nine LSBs have x 2 /N = 
2.0, indicating larger scatter than expected from observational 
errors and intrinsic scatter in the M07 relation. However, 
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since the LSBs appear to trace the M07 relation approxi- 
mately, the comparison to M31 dSphs is perhaps the most 
compelling. Masses inferr ed from the new kinematic mea- 
surem ents of M3 1 dsphs dKalirai et alJ l2009t ICollins et alJ 
2009) imply that these objects fall systematically below the 
M07 relation, with X 2 /N = 5.3. Of nine measured M31 
dSphs, all but Andl, AndVII and AndXV are outliers from 
the trend traced by spirals and MW dSphs. This behavior 
follows directly from the fact that at a given half-light ra- 
dius the velocity dispersions measured for M31 dSphs are 
systematically smaller than those of MW dSphs. The M31 
dSphs tend also to be larger than MW dSphs of similar lu- 
minosity dMcConnachie & Irwinl 12006). althoug h this trend 
appea rs to be limited to the brightest objects dKalirai et alJ 
2009). These offsets in the bulk properties of the MW and 
M31 satellites likely reflect differences in the environments 
provided by their hosts and/or the accretion histories of the 
satellites themselves. Whatever the reason, it is clear that in 
terms of the masses attributed to their dark matter halos the 
M3 1 dSphs continue to distinguish themselves from their MW 
counterparts. 



4. DISCUSSION 

M07 and W09 independently identify similarities among 
spiral and dSph dark mass profiles, respectively. M07 find that 
the rotation curves due to dark matter in spirals lie approxi- 
mately on top of each other, and W09 show that the hypoth- 
esis of a "universal" dSph mass profile is as well-supported 
by kinemat ic data as independent claims of a com mon dSph 
mass scale (Ma teo et alll993b fStrigaii et alJ2008L see below). 
Upon joining data sets, we find that spirals and MW dSphs — 
excluding outliers Her and Sgr and taking published velocity 
dispersions of the faintest satellites at face value — trace what 
appears to be a single relation under which dark mass scales 
with radius. This relation can be stated in terms of a rotation 
curve as in M07 (Equation [U, mass profile (Equation HJ, or 
constant acceleration (Equation [5]). Before discussing impli- 
cations of the new data available for M31 dSphs, we first con- 
sider the apparent scaling of MW dSph and spiral dark matter 
halos in the context of other reported scaling relationships in- 
volving these objects. 

If the M07 relation encodes a fundamental property of 
dark matter halos, then along with the appearance of a uni- 
versal dSph mass profile (W09), it contains and generalizes 
two additional resul ts reported in the recent literature. First, 
Striga riet al.l (J2008) model dSph dark matter halos and find 
that if they evaluate mass profiles at a fixed radius of 300 
pc (requiring extrapolation of the mass profiles of the small- 
est dSphs to radii of ~ 10r ha if), MW dSphs all have M 30 o ~ 
1O 7 M . If we simply evaluate the M07 mass profile at 300 pc, 
we obtain M^pp = 1. 8*}'? x 1O 7 M , instantly reproducing the 
Striga rTet al.l ([2008) result. Here it is important to note that 
the M07 relation does not imply that all dSph dark matter ha- 
los actually extend to 300 pc, or to any other radius that may 
be larger t han that of the observed tracer population. 

S econd. iKormendv & Freemanl (120041) . ISpano et al l (120081) 
and iDonato et al.l (12009) find that if they adopt cored dark 
matter halos with core radius ro and central density pp, spi- 
ral rotation curves imply a smal l range in central halo sur- 
face density. IDonato et alJ (120091) extend this analysis to MW 



dSphs and show that over 14 magnitudes in luminosity, cored 
dark matter halo models imply central dark matter surface 
densities of r^po = WO^qMq pc~ 2 . The apparent universality 
of ropo would then imply universality of the acceleration gen- 
erated by the dark matter at rp, with gpM( ro) ~ 0.5G7rropo = 
3.2+};| x 10" 9 cm s -2 dGentile et all 120091) . This is the same 
acceleration implied by the M07 relation for dark matter ha- 
los at all radii (Equation |5). Thus the M07 relation unites 
under a single scaling relation the appearances of a "univer- 
sal ma ss profile" (W09) a nd "common mass scale" for MW 
dSphs (Strigari et al. 2008]), and a "universal" central surface 
densi ty of dark matter halos ((Donato et al. 2009; G entile et"aT1 
2009). 

The new kinematic data for M3 1 dSphs indicate that these 
objects have masses systematically smaller than their MW 
counterparts, placing them below the extrapolation of the M07 
relation. In order to bring the M3 1 dSphs onto the M07 rela- 
tion, either their velocity dispersions would need to be revised 
upward by a factor of ~ 2.5 on average or their half-light radii 
(perhaps involving revision of distances) would need to be re- 
vised downward to ~ 0.3 times the current estimates. Absent 
a compelling reason to suppose that the present M3 1 data are 
grossly affected by systematic errors, it seems that the for- 
mation and/or evolution of the MW and M31 systems have 
differed sufficiently to produce measurable differences in the 
inferred properties of the dark matter halos of their dSph satel- 
lites. For example, Penarrubia et al.l (120101) use simulations to 
demonstrate that tidal interactions with the baryonic disk of a 
host galaxy reduce the masses of dSph satellites at all radii. 
These simulations reproduce the observed offset in mass be- 
tween MW and M3 1 dSphs by invoking a more-massive disk 
for M3 1 . Note that since the degree of mass loss depends on 
the parent's disk mass, this scenario requires no fine tuning of 
the orbital distributions of MW and M31 satellites. 

Left unexplained is the puzzling circumstance that the dSph 
satellites of one but not both the MW and M31 appear to 
follow the same dark matter halo scaling relation as spiral 
galaxies. Here we have emphasized the consistency of MW 
dSph data with extrapolation of the M07 relation in order 
to generalize recent and independent claims of uniformity 
among these galaxi es (Strigari et al. 2008; D onato et al.l2009t 
iGentile et alJ12009l W09). To the extent that the new data for 
M3 1 dSphs undermine the extrapolation of the M07 relation 
to dSphs, they also undermine each of these previous claims 
of uniformity. On the other hand, given the susceptibility of 
the smallest dark matter halos to evolution driven by environ- 
ment, it is feasible that one of the Local Group satellite pop- 
ulations has evolved significantly more than the other, alter- 
ing what may have been more similar conditions initially. In 
any case, it is likely that the emerging contrast between MW 
and M3 1 dSphs relates important details of the processes that 
shaped these two populations. 

This work follows directly from conversations that took 
place at the workshop Extreme Star Formation in Dwarf 
Galaxies (Ann Arbor, MI, July 2009). We are grateful to O. 
Gnedin for organizing the workshop, and to J. Penarrubia, S. 
Koposov and J. Wolf for helpful discussions. MGW acknowl- 
edges support from the STFC-funded Galaxy Formation and 
Evolution program at the Institute of Astronomy, Cambridge. 
SSM acknowledges support from NSF grant AST-0908370. 
MM and EWO acknowledge support from NSF grants AST- 
0808043 and AST-0807498, respectively. 
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